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Abstract 

Using an elementary approach, we establish a new maximum principle for the 
diffusive Lotka-Volterra system of two competing species, which involves pointwise 
estimate of an elliptic equation consisting of the second derivative of one function, 
the first derivative of another function, and a quadratic nonlinearity. This maximum 
principle gives a priori estimates for the total mass of the two species. Moreover, 
applying it to the system of three competing species leads to a nonexistence theorem 
of traveling wave solutions. 


1 Introduction 

In this paper we study the following diffusive Lotka-Volterra system of two competing 
species: 


u t = d\ Uyy + u(a 1 - c n u - c 12 v), 

yem, 

t > 0, 

Vt = d 2 Vyy + v (cr 2 c 2 i u c 22 v), 

y e M, 

t > 0, 


which is a system frequently used to model competitive behaviour between two distinct 
species. Here u(y,t ) and v(y,t ) stand for the density of the two species u and v, respec¬ 
tively; di, <7i , ca (i = 1,2), and c t j (■ i,j = 1,2 with i ^ j ) are the respective diffusion 
rates, intrinsic growth rates, intra-specific competition rates, and inter-specific compe¬ 
tition rates, all of which are assumed to be positive. The problem as to which species 
will survive in a competitive system is of importance in ecology. In order to tackle this 
problem, we consider traveling wave solutions, which are solutions of the form 

(u(y,t),v(y,t)) = (u(x),v(x)), x = y — 91, (1.2) 
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where 9 is the propagation speed of the traveling wave. In general, the sign of 9 indicates 
which species is stronger and can survive. 


We note that by using a suitable scaling, the two-species system (1.1) can be rewritten 


as 


Ut = Uyy + u(l — u — aiv), 

y eM, 

t > 0, 

v t = d v yy + k v (1 — a 2 u — v), 

y e M, 

t > 0, 


(1.3) 


where d, k, «i and a 2 are positive parameters. It is readily seen that in general, (|1.3) has 


four equilibria: ex = (0,0), e 2 = (1,0), e 3 = (0,1) and e 4 = (u*,v*), where (u*,v*) = 
( i-a°a 2 , 1 ^ 0102 ) * s ^ ie intersection of the two straight lines 1— u— cq v = 0 and 1—a 2 u—v = 
0, whenever it exists. We note that u*, v* > 0 if and only if ai, a 2 < 1 or aq, a 2 > 1. When 


the domain is bounded, the asymptotic behavior of solutions (u(y, t),v(y , t)) for (1.3) with 
initial conditions u(y, 0), v(y, 0) > 0 can be classified into four cases, as described in: 


Proposition A ([3]). Let (u(y,t),v(y,t)) be the solution of (1.3) with the entire space M. 
replaced by a bounded domain in M under the zero Neumann boundary conditions. Then 
for initial conditions u(x, 0),u(x, 0) > 0, we have 


(i) a\ < 1 < a 2 


lim (u(y,t),v(y,t)) = (1,0); 
£—>-00 


(ii) a 2 < 1 < ai 
(m) ai > 1, a 2 > 1 
(iv) ai < 1, a 2 < 1 


=> }™(u(y,t),v(y,t)) = (o,i); 

£—>■00 

(1, 0) and (0,1) are locally stable equilibria; 
=> lim (u(y,t),v(y,t)) = (■ u*,v*). 

£—>■00 


In this paper, we consider the following traveling wave problems, which are obtained 
by substituting (1.2) into (1.3) and into fll.l ) respectively, 


u xx + 9 u x + u (1 — u — Gp v) — 0, x G M, 
d v xx + 9 v x + k v (1 — a 2 u — v) = 0, i6l, 
, («, «)(-oo) = e 2 , (u, v)(+oo) = e 3 , 


(1.4) 


and 

d\ u xx + 9 u x + u (oq — di u — c 32 v) — 0, i6R, 

< d 2 v xx + 9v x + v (a 2 - c 2 iu - c 22 v) = 0, iel, (1.5) 

k (w,u)(-oo) = (£,0), («,u)(+oo) = (0, g). 

We call a solution (u(x),v(x)) of 0 an (e 2 , e 3 )-wave. The typical situation 

lim (u, v)(x) — (1,0), lim (u,v)(x) — (0,1), (1.6) 

x —>■—OO £—>-oo 


i.e. u is dominant on the left region and v is dominant on the right region of M, motivates 
us to study the (e 2 , e 3 )-wave. In this situation, u will occupy the whole domain eventually 
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if 9 > 0 while v will occupy the whole domain eventually if 9 < 0. From the viewpoint of 
ecology, we can conclude that the sign of 9 determines which species is stronger, i.e. u is 
stronger if 9 > 0 and v is stronger if 9 < 0. 

Much attention has been paid to the (e 2 , e 3 )-wave. For cases (?) or (in) in Proposition 
A, Kan-on (0,0), Fei and Carr (0), Leung, Hou and Li (|9j), and Leung and Feng 
m ) established the existence of (e 2 , e 3 )-waves employing different approaches. Under 
certain assumptions on the parameters, Mimura and Rodrigo (in mi) constructed exact 
( e 2 , e 3 )-waves by applying a judicious ansatz for solutions. By applying the hyperbolic 
tangent method, Hung([5]) found exact (e 2 , e 3 )-waves under certain assumptions on the 
parameters. All the exact (e 2 , e 3 )-waves proposed by Mimura and Rodrigo ([TO], ITT] ) 
and Hung ([5]) are represented in terms of polynomials in hyperbolic tangent functions. 
Throughout this paper, we restrict our attention to the bistable case , i.e. case (in) 
a\ > 1, a 2 > 1 in Proposition A. 

To understand the ecological capacity of the inhabitant of the two competing species, 
the investigation of the total mass or the total density of the two species u and v is 
essential since the inhabitant is resource-limited. This gives rises to the problem as to 
the estimate of u + v in (1.5). In [T|, upper and lower bounds oiu + v are given when the 


two diffusion rates d\ and d 2 are equal. However, the approach employed in [1J to obtain 
upper or lower bounds for u + v cannot be applied to the case where the diffusion rates 
d\ and d 2 are not equal. 

The above discussion raises the following questions: 


Ql: In (L.5U when d\ d 2 , can upper and lower bounds of u + v be obtained? 


As for the answer to Ql, it seems as far as we know, not available in the literature. To 
give an affirmative answer to this question, we develop a new but elementary approach. 
In fact, employing this approach leads to an affirmative answer to the following more 
general question: 


Q2: In R.5|, when d\ ^ d 2 , can upper and lower bounds of T\ u + t 2 v, where T\, t 2 > 0 


are arbitrary constants, be given? 

Since the physical units of u and v may not be identical, it makes sense to consider 
T\ u + t 2 v for the total mass in general. Although we can estimate T\ u + r 2 v via 


mm 


[ti, t 2 ] (u + v) <T 1 U + t 2 v < max[ri, r 2 ] (u + v) 


(1.7) 


once u + v is measured, more information will be wasted in this manner of estimation 
as the difference of Ti and r 2 becomes larger. Consequently, an approach which can 
accommodate to various T\ and r 2 will be of great interest. 


By adding the two equations in (1.5), we obtain an equation involving p(x) = au+/3v 
and q(x) = d±au + d 2 /3v 

0 = a [d x u xx + 9 u x + u (cti - c n u - c 12 v)) + /3 (d 2 v xx + 9v x + v (a 2 - c 21 u - c 22 v)) 

= q"(x) + 6p'(x) + au(ai — c n u- c 12 v) + /3 v (a 2 - c 21 u - c 22 v). 

( 1 . 8 ) 

The case where d\ = d 2 or p(x) is a constant multiple of q(x ) has been considered in 
[1]. Obviously, difficulties arise, and the approach used in |1] cannot be applied when 
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d\ 7 ^ d 2 since p{x) no longer can be written as a constant multiple of q(x). The approach 
proposed here can be employed to give estimates of q{x) even distinct p(x) and q(x) (i.e. 


d\ 7 ^ d 2 ) are involved in the scalar equation (3.1) 


To simplify the problem, we consider (1.4) first and present the results for (1.5) in 


Section [ 6 } One of the main results in this paper is the following maximum principle for 
Lotka-Volterra systems of two strongly competing species. 

Theorem 1.1 (Maximum Principle for q(x)). Suppose that aq > l,a 2 > 1 and 


(u(x),v(x)) is a nonnegative solution to (1.4). Then 
a f3 


mm 


mm 


[ 1 , a ] < q(x) < max 


a 


LX 3 


- &2 d CL\ - 

where q(x) = au(x) + d/3v(x) and a, (3 are arbitrary positive constants. 


max[l,d 2 ] for i 6 l, (1.9) 


In particular, we notice that the estimate of q in Theorem |1.1| does not depend on the 
propagating speed 6 and the constant k. 


The maximum principle in Theorem |1.1| can be generalized to hold true for a wider 
class of autonomous elliptic systems: 

" d\ u xx + 6 u x + u m f(u, v) = 0 , i 6 R, 

d 2 v xx + 9 v x + v n g(u, v) — 0 , ieR, ( 1 . 10 ) 

k (u,t;)(-oo) = e_, (m,u)(+oo) = e+, 

where m > 0 , n > 0 and 

e_,e+ G Cf y g = {(«, u) | u m f(u,v ) = 0 ,v n g(u,v) = 0 ,u,v > 0 }. 

We assume that f(u,v) E C°’ T (M + x M + ) and g(u,v) G C 0 , ' r (M + x M + ) for some r > 0, 
and the following property holds: 

[A] There exist u > u > 0 and v > v > 0 such that 

f(u,v ) < 0 and g(u,v) <0 if (u, v) G TZ — {(u,v) | | + | > 1, u,v > 0}; 
f(u,v ) > 0 and g(u,v ) >0 if (u,v) G TZ — {(u, v) | ^ ^ < 1, u,v > 0}. 

We have the following theorem. 

Theorem 1.2 (Generalized Maximum Principle). Assume that [A] holds. If a > 0, 
b > 0 , and (u(x),v(x)) is a nonnegative solution to ( 1 . 10 ) with e_ 7 ^ ( 0 , 0 ) and e + 7 ^ 
( 0 , 0 ), then 

( , \ nrin {d 1 ,d 2 ) . . . . , _ _x max(d 1 ,d 2 ) 

l (au, bv) -——— < au(x) + bv[x) < max ati, bv) -——— (1.11) 

v ' max(ai,a 2 ) mm(a 1 ,a 2 ) 


min 


for x G 
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Using the properties of the nonlinear terms of (1.4) more delicately, one can obtain 


better but complicated estimates for u + v. In the following, we just state an improved 
result for d — k — 1 since the form of the lower bound obtained is simple in this case. 
More general results are described in Section |4j 

Theorem 1.3. Suppose d — k — 1, aq > 1, a 2 > 1, and (u(x),v(x)) is a nonnegative 


solution to (1.4). Then for x G 


a i + 02 4 " 2 


< u(x) + v(x) < 1 . 


( 1 . 12 ) 


It is easy to see that the lower bound for u+v obtained by Theorem |l.l| is min[l/ai, l/a 2 ], 
which is smaller than or equal to +9 and is less sharp when ai,a 2 > 1. Note that 


the lower bound in ( 1 . 12 ) approaches 1 as (ai,a 2 ) approaches ( 1 , 1 ). 

As an application of Theorem we establish nonexistence of traveling waves solu¬ 
tions for the Lotka-Volterra system of three competing species, i.e.nonexistence of trav¬ 
eling solutions of 

" di u xx + 9 u x + u ( G\ — cn u — c 12 v — ci 3 w) — 0, igK, 

d 2 v xx + 6v x + v (a 2 - c 2i u - c 22 v - c 23 w ) = 0, iGK, (1.13) 

,d 3 w xx + 6w x +w (o r 3 - c 3 i u - c 32 v - c 33 w) = 0, iGK, 

where u(x,t), v(x,t) and w(x,t ) represent the density of the three species u, v and 
w respectively; cr„ Ca (i = 1,2,3), and c l3 (■ i,j = 1,2,3 ,i j ) are the diffusion 
rates, the intrinsic growth rates, the intra-specific competition rates, and the inter-specific 
competition rates, respectively. These constants are all assumed to be positive. 


For (1.13), existence of solutions with profiles of one-hump waves supplemented with 


the boundary conditions 


(u, v , w)(- oo) = ( —, 0, o) , (u, v, w)(oo) = (o, —, o) (1.14) 

is investigated in [2j. Here a one-hump wave is a traveling wave which consists of a 
forward front v, a backward front u, and a pulse w in the middle. By finding exact 
solutions and using the numerical tracking method AUTO, the existence of one-hump 


waves for (1.13),(1.14) is established under certain assumptions on the parameters ([2]). 


On the other hand, nonexistence of solutions for (1.13) and (1.14) is studied in [lj 
when the diffusion rates d\, d 2 , and d 3 are assumed to be identical. In pL|, a subtle 
structure of the competing system, which heavily relies on equal diffusivity, is employed. 
With the aid of Theorem 0 (or the extended version Theorem |6.3| in Section [ 6 ]), we give 
a much more general nonexistence of solutions for (1.13) and (1.14) when the diffusion 
rates of the species are no longer the same. 


Theorem 1.4 (Nonexistence of 3-species wave). Let 

o '2 C 33 — <73 C 23 . Assume that the following hypotheses hold: 


— oh C 33 — o 3 Ci 3 and <p 2 — 


[Hll 


> 0 ; 


5 




















[H2] C21 01 > Cn 02, C12 02 > C22 01,' 


[H3l min 


C31 02 C 32 ' 


C21 C ?2 ’ C12 <h 


mm 


[df, C?l] > CX3C33. 


Then (1.13) and (1.14) has no positive solution (u(x),v(x),w(x)). 


Biological interpretation: Due to [H2], u and v are strongly competing in (5.3) 
(see Section [5]). However, we can find parameters such that [BiS] (see the Appendix in 
Section [ 6 ]) which is slightly different from [H2] holds as well, i.e. u and v are also strongly 
competing in ( |1.13 ) as w is absent. Moreover, it is easy to see that [H3] clearly holds 
if (J 3 is sufficiently small when other parameters are fixed. In conclusion, Theorem |1.4 


asserts that, under certain conditions on the parameters, the three species u, v and w in 


the ecological system modeled by (1.13) and (1.14) cannot coexist if the intrinsic growth 


rate 03 of w is sufficiently small when strong competition between u and v occurs. 

The remainder of this paper is organized as follows. Section [2] is devoted to the proof 
of Theorem |1.1[ Then we generalize Theorem |1.1| in Section [3j By using the tangent 


line to the quadratic curve au(l — u — aiv) + (3 kv (1 — a 2 u — v) =0, it is shown in 
Section [4] that, under a certain condition on the parameters, a stronger lower bound than 
the one given in Proposition 252 and Theorem 0 can be obtained. Also, the proof of 
is presented in Section |4j As an application of Theorem 1.1, 


Theorem 1.3 


we establish 


Theorem 1.4|in Section [5} Finally, we conclude the paper with corresponding results for 


(1.5) in the Appendix (Section [b]). 


2 Proof of Theorem 11.11 


In this section p(x) = au(x) + (3v(x) and q{x) = au(x ) + d/3v(x), where a and (3 are 
arbitrary positive constants. We begin with a useful lemma. 


Lemma 2.1. Under the bistable condition a\ > 1 and ci 2 > 1, the quadratic curve 
au (1 — u — ai v) + (3 kv (1 — a 2 u — v) — 0 is a hyperbola for a > 0 and /3 > 0. 

Proof. The discriminant of the quadratic curve au (l — u — a\ v) + (3 kv (l — a- 2 U — v) = 0 is 
(a ai+/3 k a 2 ) 2 — 4 a (3 k. Since ai, a 2 > 1, we have ( a a i + 0 k CI 2) 2 > 4 a /3 ka\a 2 >4 a fdk. 
The positivity of the discriminant gives the desired result. □ 

The lemma indicates that the quadratic curve 

F(u, v ) := a u (1 — u — a\ v) + fd k v (1 — a 2 u — v) = 0 (2.1) 


cannot either be an ellipse or a parabola under the bistable condition ai, 02 > 1 . 

In Propositions |2.2 and A3 below, we give a lower bound and an upper bound for 
q(x), respectively. Combining the results in Propositions 2.2 and 2.3 
obtain Theorem 11.11 


we immediately 


6 






































Proposition 2.2 (Lower bound for q = q(x)). Let a\ > 1 and <22 > 1. Suppose that 
(u(x),v(x)) is C 2 , nonnegative, and satisfies the following differential inequalities and 
asymptotic behaviour: 

u xx + 9 u x + u (1 — u — ai v) < 0, 16I, 

< dv xx + 9 v x + kv (1 — 02 u — v) < 0, igK, (2.2) 

, (u,v) (-00) = e 2 , (u,v)(+oo) = e 3 . 


Then we have for 16 M, 


q(x ) > min 


a ft 
02 d a 1 


min[l, d 2 ]. 


(2.3) 


Proof. Let 7Z = {(w.u) 11 — u — a\V > 0,1 — a 2 u — v > 0, u > 0, v > 0}. First we construct 
an appropriate N-barrier consisting of three lines au + d/3v = A 2 , au + (3v = r) and 
au + df3v = Ai, and chose Ai, A 2 and rj as large as possible such that Q\ t C P v C Q\ 2 C 
71, where Q\ = {(u, v) \ a u + d /3v < A, u > 0, v > 0} and P v = {(u,v) \au + /3v < q,u > 
0, v > 0}. Then we show that Ai can be taken to equal the value on the right hand side 
of (2.3) and q{x) > Ai can be verified via the structure of the N-barrier. 

Now we illustrate how to construct the N-barrier in detail. For the case of d > 1 and 
f3 02 d > a a\, the N-barrier is constructed in the following three steps (see Figure 2,1(a)): 


(1) The construction of the upper pink line: we draw on the uv-plane the upper pink 

line au + d/ 3 v = X2 which passes through (-L, 0). This gives A2 = and hence the 
upper pink line is represented by the equation au + d/3v = ff. The u-coordinate 
of the u-intercept of a u + d/3 v = ^ is which is less than or equal to ^ by 

the assumption fa 2 d > aa\. This means that the u-coordinate of the u-intercept 
of au + d f3v = — is below the u-coordinate of u-intercept of the 1 — u — a x v = 0. 

(2) The construction of the yellow line: we let the the yellow line au P fdv = q start 
from (0, jffpd)- This leads to rj = and hence the yellow line is represented by 
the equation a u + f3 v = The u-coordinate of the u-intercept of a u + f3 v = 

is —bj, which is less than or equal to — by the assumption d > 1. This means that 
the u-coordinate of the u- intercept of au + f3v = ^ is less than or equal to the 
w-coordinate of w-intercept of au + d/3v = 

(3) The construction of the lower pink line: we draw the lower pink line au+d f3v = Ai 
passing through (^,0). This gives Ai = 


There are three other cases, each of which can be treated in a similar manner for the 
construction of the corresponding N-barrier (see Figures 2.1(b), 2 jl(c) and 2. 1 (d)| ) . More 
precisely, we have the following four cases and for each case, we take different Ai, A 2 and 
rj, and show that q(x) has the lower bound Ai for iGl: 


• If d > 1, 

Ct Ct Ct 

(i) when (3 a 2 d > aai, we take (Ai, A 2 , ^ 7 ) := (-—,--); 

d2 d CL 2 CL2 d 
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/. n . . ,. . . , /3 (3 d /3. 

(11 ) when p 0,2 d < a a±, we take (Ai, A2, q) := (—, —, —). 

a\ a\ a\ 

If d < 1, 

ad a a 

[in) when p 0,2 d > a a±, (Ai, A2, q) := (—, —, —); 

02 02 o2 

s , „ , , s fid 2 fid fid 

\iv) when p a 2 d < aa i, (A 1? A 2 , 77 ) := (-, —, —J. 

nx tii tlx 


We note that case (z) corresponds to Figure 2 1(a) in which the N-barrier has been 
constructed in the above three steps. The other cases (zz), (in), and (iv) correspond 
to Figures 2.1(b), 2 |1 (c) and 2.1(d), respectively. 

We first observe that the property q(x) > Ai in the four cases can be reduced to the 
following two cases: 

• for /3 02 d > a a\, q(x) > ^ min[d, 1/d] for all 16 R; 

• for j3 02 d < a ai, q(x) > — min[l, d 2 ] for all i 6 l 


Combining the two cases above leads to q(x) > min[^, p] min[l, d 2 ] for all x E M, which 
is the desired result. 

Now we show q(x) > Ai in (z) ~ (iv). The two inequalities in (2.2) and (2.1) give 

q”(x) + 9p'(x) + F(u(x), v(x)) <0. (2.4) 

For d > 1, we first prove (i) by contradiction. Suppose that, contrary to our claim, 
there exists zeR such that q(z) < Ai. Since u,v G C 2 (M), by (u,v)(— 00 ) = (1,0) and 
(m,z:)(+oo) = (0,1), we may assume min xe K q(x) = q(z). We denote respectively by Z2 
and Z\ the first points at which the solution (u(x),v(x)) intersects the line au + d f3v = X 2 
in the on-plane when x moves from z towards 00 and —00 (as shown in Figure 2.1(a)). 
For the case where 6 < 0, we integrate (2.4) with respect to x from z 1 to z and obtain 

q(z) — q'(zi) + 0 (p(z) — p(z{]) + f F(u(x),v(x)) dx < 0. (2.5) 

J Zi 

On the other hand we have: 


• since min xe R q(x) = q(z), q'(z) = au'(z) + d/3v'(z) = 0 ; 

• q(z\) = A 2 follows from the fact that Z\ is on the line au + d/3v = A 2 . Since z\ 
is the first point for q(x) taking the value A 2 when x moves from z to — 00 , we 
conclude that q(z\ + 5) < A 2 for z — Z\ > 8 > 0 and q'(z{) < 0; 

• p(z) < r] since z is below the line au + (3v = iy, p(z{) > q since z is above the line 
au + j3v — rj; 


it is readily seen that the quadratic curve F(u,v) = 0 passes through the points 
(0,0), (1,0), (0,1), and (u*,v*) in the tm-plane. Let A + = {(u,u) | F(u,v) > 0 ,u> 
0,n > 0}. By Lemma 2.1 and the property that F(u,v) < 0 for large u and v, it 
follows that A + is the region bounded by a hyperbola, w-axis and w-axis. Moreover, 
{(m(x), v(x)) I z\ < x < zj C TZ C A + . Therefore we have f~ F(u(x),v(x)) dx > 0. 
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Summarizing the above arguments, we obtain 


q'(z) ~ q'(zf) + 9 (p(z) — p(zi)) + / F(u(x),v(x)) dx > 0, 

J Z\ 


( 2 . 6 ) 


which contradicts (2.5). Therefore when 9 < 0, q(x ) > Ai for x G M. For the case where 


9 > 0, integrating (2.4) with respect to x from z to z 2 yields 


"Z2 


q'(z 2 ) - q'(z) + 9 (p(z 2 ) — p(z)) T / F(u(x),v(x)) dx < 0. (2.7) 

J Z 

In a similar manner, it can be shown that q'(z 2 ) > 0, q'(z) = 0, p(z 2 ) > V, p(z) < V, and 


Jf 2 F(u(x),v(x)) dx > 0 . These together contradict (2.7). Consequently, («) is proved for 


d > 1. For d — 1, we have q = p and (2.4) becomes 


p"(x) + 9p(x) + F(u(x),v(x)) < 0, 16 


( 2 . 8 ) 


Moreover, when d — 1 we take Ai = A 2 = rj = —, i.e. the three lines au + d/3v = Ai, 
au + d/3v = A 2 , and au + j3v = q coincide. Analogously to the case of d > 1, we 
assume that there exists z G M such that p(z) < Ai and min xS Rp(x) = p(z). Due to 
min a , e Mp(x) = p(z), we have p'(z) = 0 and p"(z) > 0. Since {u(z),v(z)) is in the interior 
of TZ, which is contained in the interior of A + , we have F(u(z),v(z)) > 0. These together 
give p"(z) + 9p\z ) + F(u(£),v(z)) > 0, which contradicts (2.8). Thus, p(x) > Ai for all 
x E I when d—l. As a result, the proof of {%) is completed. 

The proofs for cases (ii), (in), and (■ iv ) are similar (see Figures 2.1(b), 2.1(c), and 
2,l(d)[). This completes the proof of Proposition 2.2 


□ 


Proposition 2.3 (Upper bound for q = q(x)). Assume that a\ > 1, a 2 > 1, and that 
(u(x),v(x)) is C 2 , nonnegative, and satisfies the following differential inequalities: 

u xx + 9u x + u (1 — u — a\ v) > 0, x G M, 

< dv xx + 9v x + kv (1 — a 2 u — v) > 0, x G M, (2.9) 

, (u, v)(-oo) = e 2 , ( u , v) (Too) = e 3 . 

Then for x G M, we have 


q(x) < max 


,P 


max[l, d 2 


( 2 . 10 ) 


Proof. As in the proof of Proposition 2.2 there are also four cases and for each case, 


we can construct the N-barrier as shown in Figures 2.2(a), 2.2(b). 2 |2(c) and 2.2(d) and 
prove that q(x) < Ai for x G M: 


• If d > 1, 

(i) when /3 d > a, we take (Ai, A 2 , rj) := (j3 d 2 , /3 d, (3 d); 

(ii) when (3 d < a, (Ai, A 2 , rj) := (a d, a, a). 
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v a- "ta s 




(a) (b) 




(c) (d) 


Figure 2.1: Red line: 1 — u — a\ v = 0; blue line: 1 — 02 u — v = 0; green curve: au(l — u — 
01 v) + /3 k v (1 — 02 u — v) = 0; magenta line (above): au + d f3v = A 2 ; magenta line (below): 
au + d/3v = Ai; yellow line: au + /3v = 77 ; dashed curve: (u(x),v(x)). |(a)| a\ = 2, 02 = 3, 


= 17, /3 = 18, and d = 2 give Ai = A 2 = and rj = (b) ai = 2, 02 = 3, a = 17, 

= 5, and d = 2 give Ai = §, A 2 = 5, and r/ = |. (c) ai = 2, 02 = 3, q = 17, /3 = 18, and 


= | give Ai = 


34 
9 ’ 


A 2 = ^, and 77 = 


17 

3 


= f, A 2 = and 77 = 


(d) ai = 2, 02 = 3, a = 17, /3 = 18, and d = \ give 
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• If d < 1, 


(Hi) when (3 d > a, (Ai, A 2 , rj) := (/ 3 , /3 d, (3)] 
(iv) when (3 d < a, (Ai, A 2 , rj) := (^, a, ^). 


We note that cases (i), (ii), (in), and (iv) corresponds to Figures 2.2(a), 2.2(b), 2.2(c) 
and 2.2(d), respectively. Combining the four cases above, it follows that 


for (3 d > a, q(x) < (3 max(l, d 2 ) for all x E 


• for (3 d < a, q(x) < a max(d, 1/d) for all 16 R, 


which implies q(x) < max^,^] max[l,<i 2 ] for all 1 6 1 
similar to that of Proposition 222 and is hence omitted. 


The rest part of the proof is 

□ 


3 General maximum principle 


In this section, we prove Theorem 1.2 which generalizes the maximum principle in 


Theorem |1.1| to elliptic systems with a wider class of nonlinear terms. Recall that 

*,5)1 


n = {(w,r07 5 + I > !, u,v > 0} and K = {(u,v ) \ \ + £ < 1, u,v > 0}. 

As in Section [2j adding the two equations in (1.10) leads to an equation involving 
p(x) = a u(x) + [3 v(x) and q(x) = d\a u(x) + d 2 (3v(x), i.e. 

0 = ol ( diu xx + 9u x + u m f(u,v)) + f3 (d 2 v xx + 9v x + v n g(u,v)) 

= q”(x) + 9p'(x) + au m f(u, v) + f3 v n g(u, v) 

= q"(x) + 9p'(x) + F(u,v), (3.1) 

where a, (3 > 0 are arbitrary constants and F(u,v) = au m f(u,v) + f3v n g(u,v). By 
assumption [A], it readily follows that F(u,v) > 0 on 1Z and F(u,v) < 0 on TZ. In 
when e_ = (0,0) or e + = (0,0), the lower bound estimate no longer holds 


1.2 


Theorem 

but the upper bound estimate is still valid. In the following, we state a theorem, which 


is slightly more general than Theorem |1.2[ to include the upper bound estimate when 
e_ = ( 0 , 0 ) or e + = ( 0 , 0 ). 

Theorem 3.1. Assume that [A] holds. If a > 0, b > 0, and (u(x),v(x)) is a nonnegative 


solution to (1.10), then 


where 


and 


A < au(x) + bv(x) < A, 16 I, 


T { ~ U~\ lllaX 

A = max ( au,bv )- 

v ' mini 

min(o 


1 , bv) 
A = min (au,bv) 


max( 


(3.2) 

(3.3) 

(3.4) 
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OX) 0.5 L0 1.5 2.0 


u 

(c) 



(d) 


Figure 2.2: Red line: 1 — u — au v = 0; blue line: 1 — 02 u — v = 0; green curve: au(l — u — 
01 v) + /3 k v (1 — 02 u — v) = 0; magenta line (below): au + d (3v = A 2 ; magenta line (above): 
au + d/3v = Ai; yellow line: au + (3v = rj; dashed curve: (u(x),v(x)). |(a)| a\ = 2, 02 = 3, 
a = 17, j3 = 18, and d = 2 give Ai = 72, A 2 = 36, and r\ = 36. |(b)| a\ = 2, 02 = 3, a = 17, 
[3 = 5, and d = 2 give Ai = 34, A 2 = 17, and 7 / = 17. (c) au = 2, 02 = 3, a = 17, /3 = 33, and 
d = | give Ai = 33, A 2 = 22, and r/ = 33. (d) au = 2, au = 3, a = 17, /3 = 18, and d = | give 
Ai = 34, A 2 = 17, and rj = 34. 
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with x defined by 


X = 


0, if e+ = (0,0) or e_ = (0,0), 
1, otherwise. 


(3.5) 


Proof. Let a = a/d\ and (3 = b/d 2 in (3.1). Then q(x) = au(x) + bv(x) in (3.1). We 
employ the N-barrier method developed in Section[2]to show (3.2), which implies (1.11). 

First, we assume e + 7 ^ (0,0) and e_ 7 ^ (0,0). To construct an appropriate N-barrier 
for the lower bound estimate, we consider Q\ = {(u, v) \d\ a u + d 2 (3 v < A, u > 0, v > 0} 
and P v = {(«, v) | au + (3v < rj, u > 0, v > 0}, and chose Ai, X 2 and rj as large as possible 


such that Q Xl C P^ C Q\ 2 
by 


C TZ. By direct computation, Ai, A 2 and rj can be determined 


A 2 = min (diau, d 2 ay) = min (an, by '), 

■ f 1 

77 = mm ( —, —) A 2 , 
a 1 a 2 

Ai = min (d 1 ,d 2 ) rj = min (^-, ^) A 2 . 

Cfc2 d\ 


(3.6) 

(3.7) 

(3.8) 


Since F(u,v) > 0 on TZ, we can employ (3.1) and follow the arguments in Section [ 2 ] to 
obtain the lower bound estimate q{x) > \\ = \. 

When the boundary conditions e + = (0,0) or e_ = (0,0), the argument in Section [2] 
can not be applied and only a trivial lower bound au(x) + bv(x) > 0 can be given. 

The proof for the upper bound of au + bv is similar. □ 


4 More delicate lower bound: tangent lines to quadratic 
curves 


In this section we provide an alternative approach to determine the line au + d(3v = X 2 


in the proof of Proposition |2.2| so that a bigger A 2 can be chosen and a stronger lower 
bound for q(x) can be given. To this end, we determine A 2 by solving 


a u + d(3v = A 2 , 


(4.1a) 


F(u,v) = 0, 
F u (u,v) F v (u,v ) 


a 


d /3 


(4.1b) 


(4.1c) 


where F u (u, v ) = (n, v) and F v (u, v ) = ^y(u, v). In (4.1), a u+d/3 v = A 2 is the tangent 


line to the quadratic curve F(u,v ) = 0 and the line a u + d(3v = X 2 is perpendicular to 


the vector < F u (u,v), F v (u,v) >. The solution (u,v, X 2 ) of (4.1) determines the point of 


tangency of the line au + d/3v = X 2 and the quadratic curve F(u, v ) = 0. 


Equation (4.1) can be solved with the aid of Mathematica. It is easy to see that the 


first and third equations in (4.1) are linear, while the second one is quadratic. To solve 
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(4.1), we begin by solving the first and third equations to obtain (u,v) = (w(A 2 ), v(A 2 )) 
as , namely 

—dA 2 (aai + a 2 (3k) + aj3d(d — k) + 2ak\ 2 . . 

U — -— 7 — -„ , N -^ -- ? y^.ZdbJ 


V = 


2a (—d (aai + a 2 (3k) + (3d 2 + ak) 

OM\ A 2 T /3 (a 2 k\ 2 T ad — 2dA 2 — a/c) 
2/3 (aaid + a 2 (3dk — (3d 2 — ak ) 


(4.2b) 


Substituting (4.2) into the second equation in (4.1) yields the following quadratic equation 
for A 2 : 


where 


d 2 = 


fh A 2 + /^i A 2 + fio — 0, 
(aai + CL 2 (3k ) 2 — 4 a(3k 


iii — 


4a(3 (—d (aai + a 2 (3k) + /3d 2 + afc) ’ 

—2a/3 ((<3 + fc) (aai + a 2 (3k ) — 2 k(a + ^d)) 
4a/3 (—d (aai + o, 2 (3k) + (3d 2 + afc) 

a/3(d — A:) 2 


do — 


It follows from (4.3) that 


4 (—d (aai + a 2 (3k) + /3d 2 + afc) 
Hi ± y/Hi ~ 4 h 0 h 2 


A 2 — 


2/i 2 


Using (4.4), the discriminant T> of (4.3) is given by 

k (—aia — a 2 (3k + a + (3k) 


V = Hi ~ 4 d2do = 


-d (aai + a 2 (3k) + (3d 2 + ak 


(4.3) 

(4.4a) 

(4.4b) 

(4.4c) 

(4.5) 

(4.6) 


To apply the approach proposed here, it is necessary that V > 0. In fact, D ^ 0 since 
k (—aia — a 2 (3k + a + (3k) = k(a( 1 — ai) + k(3(l — a 2 )) < 0. Moreover, it can be shown 
that T> > 0 if and only if 

aai + a 2 (3k — J (aai + a 2 (3k) 2 — 4a(3k aai + a 2 (3k + J (aai + 02 (3k) 2 — 4a(3k 

- < d < - 


2/3 


2/3 


(4.7) 

Under the condition (4.7), Ho < 0, Hi > 0, and /i 2 < 0 and hence the two roots A 2 given 
by (4.5) are both positive. However, when 

~Hi + \/Hi ~ 4 ^ 


Ao — 


2^2 


(4.8) 


it turns out that one of u,v given by (4.2) is negative. We remark that this fact can also 
be easily seen from a property of the hyperbola F(u,v) = 0. That is, for a given slope, 
there exist two tangent lines to the hyperbola F(u, v) = 0: one has the intersection point 
in the first quadrant, while the other has the intersection point in the second or fourth 
quadrant. Therefore, we have 


A 2 — 


~Hi — \J h\ — 4/j 0 /J2 
2/12 


(4.9) 
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and the intersection point (u,v) can be expressed in terms of A 2 using (4.2). Now we are 
in a position to prove: 


Proposition 4.1 (Stronger lower bound for q = q(x)). Assume that a\ > 1, 02 > 1 

and the condition (4.7) holds. Let X 2 be given by (4.9) and let (u(x),v(x)) be a pair of 
nonnegative C 2 functions satisfying the differential inequalities (2.2). Then we have 

X 2 

• when d > 1, q(x) > — for all x G R; 

• when d < 1, q{x) > A 2 d for all 16 R. 

Consequently, for d > 0 we have q(x) > X 2 min[d, 1/d] for all x G R. 


Proof. The proof is similar to that of Proposition 2.2 (see Figure 4.1 and Figure 4.2). 
Here we only explain how to determine (Ai, 77 ) (see Figure 4.2). When d < 1, we choose 
(Ai, 77 ) = (A 2 d, A 2 ) via the following procedure. The line au + d/3v = X 2 intersects the 
■u-axis and u-axis at (^,0) and (0, ^), respectively. Since d < 1 , the line au + (3 v = q 

= A2. 


passes through the point ( 77 , 0 ), and q is determined by q = au + /3v 

\2 
P 


The line au + f3v = q intersects the u-axis at (0, so Ai is given by Ai = au + 


d/3v 


x = X 2 d. For the case d > 1, Ai and q can be determined in a similar 
(«,t 0 =(o,£) ' 


manner as (Ai, 77 ) = ^f). (Refer to Figure 4.1 ) 


□ 


Combining Propositions 2.3 and 4.1, we obtain 


Corollary 4.2 (Maximum principle for q(x)). Assume that a± > 1, a 2 > 1 and 

the condition ( 4.7[ ) holds. Let X 2 be given by (4.9) and let (u(x),v(x)) be a nonnegative 
solution to the following differential equations and asymptotic conditions 


' u xx + 6 u x + u (1 — u — ai v) = 0, 

igR, 


v xx + 9 v x + v (1 — a 2 u — v) = 0, 

X G R, 

(4.10) 

k (u,u)(-oo) = (1,0), (u, u)(Too 

= (o, 1 ). 


have 



X 2 . ,2n / \ f a ~ 

— mm[l,a J < q(x) < max —,p 

max[l, d 2 ]. 

(4.11) 


Then for x G 


Compared with Theorem H3 Corollary |4.2| asserts that when Q4.7[) h olds, a stronger 
lower bound for q(x) can be given in terms of A 2 , which is dehned by (4.9). I 11 particular, 


when a = /3 = d = k = 1, we obtain Theorem 1.3 


To illustrate Proposition |4.1[ we give an example. When a\ = 2, a 2 = 3, a = 17, 
/3 = 18, d — 2, and k — 1, (4.1) can be solved to give 

/3 (2349 ± 71V4611) 17 (1131 ± 4^4611) 153 (79 ± V4611) A 
' U,Vl “ V 47270 ’ 141810 ’ 1630 J’ ' * 
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which is approximately (u,v, A 2 ) = (—0.157,0.103,1.0415) or (0.455,0.168,13.789). We 
choose (u.,..A 2 ) = ggl , "‘’’S.r 11 . 153(7 7 3 f n) ) « (0.455.0.168.13.789). 


and determine 


Al = "=7 


(4.13) 


by employing Proposition |4.1 Then we are led to Figure |4~T Applying Proposition 4.1| 
again, it follows that q{x) > Ai ~ 6.895 for all (see Figure [47T| ) . This lower bound 

is much bigger compared with the one given previously in Section [2j where the lower 
bound for q(x) is A? « 2.833 (see Figure 2.1(a)[). 



Figure 4.1: Red line: 1 — u — a± v = 0; blue line: 1 — 02 u — v = 0; green curve: 
au (1 — u — a\ v) +/3 kv (1 — 02 u — v) = 0; magenta line (above): au + d/3v = A 2 ; magenta 
line (below): au + d/3v = Ay yellow line: au + f3v = 77 ; dashed curve: (u(x),v(x)). 

ai = 2, 02 = 3, a = 17, /? = 18, d = 2, k — 1, and A 2 = ~ 13.789 give 

Ai = 77 = 1 ~ 6.895 according to (14.131) . 


16 

























Figure 4.2: Red line: 1 — u — a±v = 0; blue line: 1 — a 2 u — v = 0; green curve: 
au (1 — u — a± v) +/3 kv (1 — 02 u — v) = 0; magenta line (upper): au + d/3v = X 2 ; magenta 
line (lower): au + d/3v = Ai; yellow line: au + /3v = 77 ; dashed curve: (u(x), v(x)). 

ci\ = 2, 0,2 = 3, a = 17, /3 — 18, d — k — 1, and A 2 = 51 1 133 ^^) 1 ^ ^) ~ 8.126 give 

Ai = A 2 d = XL (133 + a/16059) w 5.418 and 77 = A 2 = r,1 ( 133 +^ 5M ) ~ 8 .i 2 6 . 


5 Application to the nonexistence of three species 


travelling waves: proof of Theorem 1.4 


In this section, we prove Theorem 1.4 by contradiction. 


Proof of Theorem l.f . Suppose to the contrary that there exists a solution (u(x), v(x), w(x)) 
to (1.13),(1.14). Due to the fact that w(x) > 0 for x G M and iu(±oo) = 0, we can hnd 


xq e M such that max l£ iui(i) = w(x o) > 0, w"(x o) < 0, and w'(xq ) = 0. Since w(x) 
satisfies d% w xx + 6 w x + w{a^ — c^\u — C32 v — C33 w) = 0, we obtain 


which gives 


w 


o '3 - c 3 i u(x 0 ) - c 32 v(x 0 ) - c 33 w(x 0 ) > 0 , 


(x) < w(x 0 ) < — (<r 3 - c 3 i u(xo) - c 32 v(x 0 )) < —, 
C 33 c 33 


x G 


(5.1) 


(5.2) 
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As a consequence, we have 


di u xx + 9 u x + w(oq - ci 3 a 3 C33 1 - cn u - c 12 v) <0, x G 
d 2 v xx + 9 v x + v(a 2 - c 23 a 3 c^ 1 - c 21 u - c 22 v) < 0 , x € 


(5.3) 


Because of [HI] and [H2], we can apply Proposition 6.1 from the Appendix to 
Indeed, [HI] assures the positivity of a 1 — C 13 a 3 c^ and a 2 — c 23 cr 3 while the 


(5.3). 


bility condition [BiS] (see Appendix) for the nonlinearity in (5.3) follows from 
Consequently, we obtain a lower bound of c 3 i u{x) + c 32 v(x), i.e. 


bista- 

[H2l. 


C31 u(x) + c 32 v(x) > C33 1 min 


C3I 02 C 32 01 
C 21 d 2 C 12 d,\ 


nun 


in [cZf, , igR. (5.4) 


The condition [H3] then yields 

c 3 i u(x) + c 32 v{x) > (T3, x e 


which contradicts (5.1). This completes the proof. 


(5.5) 


□ 


6 Appendix 


After suitable scaling, system (1.4) is equivalent to (1.5). Theorem 0 establishes lower- 
upper bound estimates for (1.4). In this section, we state corresponding results for (1.5). 


Throughout this section, we shall always assume the bistable condition: 

[BiS] > £2-, ^ 

L J cn C21 ’ C22 C12 ’ 

which correponds to the condition a\ > 1 and a 2 > 1 used in previous sections. Let 


p(x) = au + /3v and q(x) = d 3 au + d 2 (3v. Then, from (1.5), it follows that p(x) and 
q(x ) satisfy 

a {di u xx + 9u x + u (oq -c u u- c 12 v)) + 0 (d 2 v xx + 9v x + v (cr 2 -c 21 u- c 22 v)) 

= q"(x) + 9p'(x) + au (oq - c u u - c 12 v ) + /3 v (cr 2 - c 21 u - c 22 v). 

( 6 . 1 ) 

We can now apply the approach proposed in Section [2] to obtain lower and upper bounds 
for q(x) in Proposition 6.1 and Proposition 6.2 respectively. 


Proposition 6.1 (Lower bound for q = q(x)). Suppose that (u(x),v(x)) is C 2 and 
nonnegative, and satisfies the differential inequalities 

di u xx + 9 u x + u (<Ti — cn u — C 12 v) <0, 16 K, 

< d 2 v xx +9v x + v (cr 2 - c 2 i u - c 22 v) < 0 , 16 I, ( 6 . 2 ) 

(u,v)(-oo) = (gL,0), (m,u)(+cx 3) = (0,g). 
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Then for x E M, we have 


q(x) > min 


a a 2 (3 a i 


C 21 d-2 c 12 d± 


mm 


in [<if, ti|] 


(6.3) 


Proposition 6.2 (Upper bound for q = q(x)). Suppose that (u(x),v(x)) is C 2 and 
nonnegative, and satisfies the following differential inequalities and asymptotic conditions 


d\ u xx + 0 u x + u (ci — cn u — C 12 v) >0, ieR, 
d 2 v xx +9v x + v ( a- 2 c 2 i u c 22 v) >0, i6l, 
(u,»)(-oo) = (gpO), (m,w)(+oo) = (0,g). 


(6.4) 


Then for x eR, we have 

q(x) < max 


a a i /3 cr 2 


c n d 2 c 22 d\ 


max 


[d\, dl]. 


(6.5) 


Theorem 6.3 (Maximum principle for q(x)). Suppose that (u(x),v(x)) is a nonneg¬ 
ative solution to the differential equations 


di u xx + 9 u x + u (a i — Cn u — ci 2 v) = 0, i6l, 
da v xx +9v x + v (cr 2 - c- 2 i u - c 22 v ) = 0, i 6K, 
(u,u)(-oo) = (g-,0), (u,u)(+oo) = (0, g|). 


( 6 . 6 ) 


Then for x E M, we have 
a a 2 f3 cri 


nnn 


C 2 i d 2 Ci 2 di 


mm 


[dl,d 2 2 \ < q(x) < 


max 


acri (3 cr 2 


Cn d 2 c 22 di 


max 


[dl,d 2 2 \. (6.7) 


Proof. Combining Proposition 6.1 and Proposition |6.2[ Theorem 6.3 is established. □ 


We note in particular that, when a = d 1 1 and (3 = d 2 1 , (|6.7[) becomes 


7-1 


mm 


-1 

1 ^ 

1 10 

3 

min 

d\ d 2 

< m(x) + n (ic) < max 

O' l cr 2 

max 


_C21 Ci 2 _ 




.di C 22. 


C?2 


( 6 . 8 ) 


The above result can be generalised by letting a = d 1 1 rq and (3 = d 2 1 r 2 . Equation (6.7) 
then leads to 


mm 


ricr 2 

dCi 

min 

di 

d 2 

. di 

Cl2 . 

d 2 

’ di_ 


< riu(x) + r 2 v(x ) 
Tier i r 2 a 2 


< max 


c 11 c 22 


max 


(l | d 2 
d 2 di 


(6.9) 
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where r 1 , r 2 > 0 are arbitrary constants. 
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